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CLAIMS! 

1. A method for processing a digital image Bi, the image Bj being a 
convolution, of an image F and a point spread function h t comprising removing 
noise from the image Bi so as to produce an image B' of reduced noise, and 

5 calculating F based upon B '. 

2. The method of Claim 1 wherein an amount of noise is calculated in a 
plurality of images B, and the image B' is selected as an image of essentially 
minimal noise among the images B. 

3. The method of Claim 1, wherein the- amount of noise in an image is 
10 calculated using an algebraic expression involving the gradient of a function P(x) 

obtained by inverse Fourier transform of B(q) I h{q) . 

4. The method of Claim 3, wherein the amount of noise N in an image B is 
calculated according to the algebraic expression N = frp(x)- VP*(x)dx t wherein 
V indicates the gradient and "*" indicates complex conjugate. 

•\ 5 5. The method according to Claim 4 wherein B\q) , the Fourier transform of 
B\ is equal to % •*-)(<?) for same integer i, where is obtained according to 

the algebraic expression B~ f + i(q) = B\(q)Q.+ g^p(qj^ q 2 )' , where s is a small 
positive number. 

6. The method according to Claim 4 wherein B\q) is obtained according to 
20 the algebraic expression B\q) = B\{qy a ^' m ^ " ,where a is a predetermined 

constant, and and, Dig) is the Fourier transform of 1/h. 

7. The method according to Claim 1 wherein calculating F involves 
calculating an inverse Fourier transform of the algebraic expression B\q)/h{q), 
wherein B'(q) is the Fourier transform of the image B' of reduced noise, and k (g) is 

25 the Fourier transform of h. 
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8. A metiiod for processing a deconvolved image B, the image B having 
been deconvolved according to a deconvolution filter D, the method comprising 
reducing correlation between the image and the deconvolution filter. 

9. The method of Claim 8 wherein an amount of correlation is calculated in a 
5 plurality of images P, and an image P' is selected among the images P as an image 

having essentially minimal correlation -with the deeonvlution filter. 

10. The method of Claim. 9 wherein the amount of correlation C in an image P 
is. calculated according to the algebraic expression C = J%jl)(£)| 2 'p(?)f wherein 
P(q) is th e Fourier transform of an image P. 

io 11. The method according to Claim 9 wherein the Fourier transform of 
P\ is equal to Pi*\{q) for same integer i, where Pi+i(q) is obtained according to 
the algebraic expression. P, + r.(g) = Pi{q)Q. + 4^Mf ^ ' where 5 * s a sma11 P ositive 
number, 

12. The method according to Claim 9 wherein P\q) is obtained according to 

15 the algebraic expression F (q) = P^q)e~^ D{q ^ ,where (3 is a predetermined constant. 

13. The method for processing a digital image Bi, the image Bi being a 
convolution of an image F and a point spread function h comprising the steps of: 

(a) removing noise from the image Bi so as to produce an image B' of 
reduced noise; 

20 (b) obtaining function P^q) according to the algebraic expression 

(c) reducing correlation between ^and l/hso as to product a function 

P' of reduced correlation; and 

(d) obtaining a rectified image F by inverse Fourier transform of P\q) . 

25 14. A method for obtaining a radius r of a point spread function h describing 
an out-of-focus distortion of a digital image B, the method comprising a step of 
calculating a gradient at a plurality of pixels in the image B. 
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15. Hie method according, to Claim 14 m which a radius r(x) is calculated at 
each of the plurality of pixels based upon the gradient 

16. The method according to Claim 15 wherein each of the plurality of pixels 
is located at an edge of the image B. 

17. The method according to Claim 15 wherein a radius r(x) is inversely 
proportional to the gradient at x. 

18. The method according to Claim 16 wherein r is obtained as the r(x) having 
an essentially maximal frequency among the calculated radii r(x). 

19. The method according to Claim 18 wherein a radius r(x) is calculated 

according to the algebraic expression '(x)-^, wherein 5 < x ) is fce absolute 
value of the gradient of B at x normalized by dividing by the height of the edge 
atx. 

20. The method according to Claim 1, furfher comprising a step of producing 
the image B» from an image Bq, where the image B 0 was obtained using a digital 
camera that applies a transformation to a light level detected at a pixel, the 
transformation having an inverse, wherein Bi is obtained from the image B 0 by 
applying to the image B 0 the inverse transformation. 

21. A program storage device readable by machine, tangibly embodying a 
program of instructions executable by the machine to perform method steps for 

j processing a digital image B h the image B t being a convolution of an image F and 
a point spread runction h, the method comprising removing noise from the image 
Bi so as to produce an image B' of reduced noise, and calculating F based upon B'. 

22. A program storage device readable by machine, tangibly embodying a 
program of instructions executable by the machine to perform method steps for 

5 processing a deconvolved image B, the image B having been deconvoluted 
according to a deconvolution filter D 5 the method comprising reducing correlation 
between the image and the deconvolution filter. 

23. The program storage device readable by machine, tangibly embodying a 
program of instructions executable by the machine to perform method steps for 



processing a digital image Bi, the image Bi being a convolution, of an image F and 
a point spread function k comprising the steps of: 

(a) removing noise from the image Bi so as to produce an image B 1 of 
reduced noise; 

(b) obtaining function P } (q) according to the algebraic expression 

P ] (q) = B'(q)/h(q) ; 

(c) reducing calculation between Ijjand T/k&o as to product a function 

P' of reduced correlation; and 

(d) obtaining a rectified image F by inverse Fourier transform of F(q) . 

24. A program storage device readable by machine, tangibly embodying a 
program of instructions executable by the machine to perform method steps for 
obtaining a radius r of a point spread function h describing an out-of-focus 
distortion of a digital image B, the meiihod comprising a step of calculating a 
gradient at a plurality of pixels in. the image B. 



